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A symmetric spin-1/2 Ising-Heisenberg diamond chain 
with the Ising four-spin interaction is exactly solved by 
means of the generalized decoration-iteration mapping 
transformation. The ground state, the magnetization pro- 
cess and thermodynamics are particularly examined for 
the case of antiferromagnetic pair interactions (Ising and 
isotropic Heisenberg ones). It is shown that an interplay 
between pair interactions, the four-spin interaction and 



the external magnetic field gives rise to several quantum 
ground states with entangled spin states in addition to 
some semi-classically ordered ones. Besides, the tem- 
perature dependence of the magnetic susceptibility mul- 
tiplied by the temperature is studied and the interesting 
triple-peak specific heat curve is also detected when con- 
sidering the zero-field region rather close to the triple 
point, where three different ground states coexist. 
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1 Introduction Spin systems with multispin ex- 
change interactions represent objects of scientific interest 
in the past few years. Research in this field leads to a deeper 
understanding of many interesting physical phenomena, 
such as the non-universal critical behaviour |Q]|2], optical 
conductivity Raman peaks H, as well as, deviations 
from the Bloch T 3 / 2 law at low temperatures ||5]|6) . More- 
over, the effect of magnetic field on the ground state of 
quantum systems with cyclic four-spin interaction has been 
recently particularly examined as well [7l8i9 |. Of course, 
the immense theoretical interest to the spin models with 
multispin interactions is not purposeless. Multispin inter- 
actions have been experimentally observed in real triangu- 
lar magnetic systems composed of 3 He atoms absorbed on 
graphite surfaces ITOl . the hydrogen bonded ferroelectrics 
PbHP0 4 and PbDP0 4 HQ, as well as, the squaric acid 
crystal (H 2 C 2 4 ) H12lll3l[T4l and some copolymers ifTBI . 
Recently, it was shown that the cyclic four-spin interac- 
tion could also explain the neutron-scattering experiments 
concerning high-T c compounds such as La 2 Cu0 4 lfl6l . 
La 6 Ca 8 Cu 24 4 i Q2EDD, and La 4 Sri Cu 24 O 4 i fl9l . 



On the theoretical side, in spite of the numerous nu- 
merical studies predicting rich phase diagrams |2() 2 1 II22I 
l23l . a number of important questions concerning the spin 
ordering and the critical behaviour of quantum Heisen- 
berg models realized by the four-spin interaction, remain 
unclear due to controversial results obtained by different 
treatments 124,25,26,271. From this point of view, the 
exactly solvable models play an important role in under- 
standing the multispin effects. Of course, the quantum 
spin models is very difficult to deal with exactly due to a 
rather cumbersome and sophisticated mathematics, which 
precludes an exact treatment of the most (even simple- 
minded) spin systems. Motivated by this fact, a special 
class of the hybrid Ising-Heisenberg models [ 28 29 301I3T1 
fMEMMMMMMMEMIEM^ has been recently 
proposed, which overcome the afore-mentioned mathe- 
matical difficulty by introducing the Ising spins at nodal 
lattice sites and the Heisenberg dimers on interstitial dec- 
orating sites of the considered lattice. These simplified 
classical-quantum spin models can be rigorously solved 
by two different exact analytical approaches; the clas- 
sical transfer-matrix technique [44,45], which is rather 
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Figure 1 A part of the spin- 1/2 Ising-Heisenberg diamond chain 
with the four-spin interaction. The empty (full) circles denote 
lattice positions of the Ising (Heisenberg) spins. The ellipse de- 
marcates spins belonging to the fcth diamond unit. 



straightforward but useable only for one-dimensional spin 
systems, and/or more universal generalized decoration- 
iteration transformation [46 47 48 49 1, which can be used 
for both the one- and two-dimensional systems. It is also 
worth mentioning that in order to investigate many inter- 
esting physical phenomena, various extensions and gener- 
alizations of the hybrid Ising-Heisenberg models may be 
done without loss of their exact solubility. For instance, 
it is possible to extend the Ising-Heisenberg models by 
including the next-nearest-neighbour exchange interaction 
between the Ising spins 1501 . the Dzyaloshinskii-Moriya 
anisotropy acting on the decorating Heisenberg spins BP . 
or to solve exactly the analogous Ising-Heisenberg mod- 
els with the Heisenberg spins S > 1/2 |29]|35]|52]|53]|53 
[55l[5M57l|58ll59ll60ll6T1l that bring a deeper insight into 
how the magnetic behaviour of the quantum spin systems 
depends on the magnitude S |35j56J[57| of the decorat- 
ing Heisenberg spins and also allow to examine the ef- 
fect of other interaction terms such as the axial zero-field 
splitting parameter J29]|52|53l|5l|55l|58l|59]|60l and/or the 
biquadratic XXZ interaction 160 )611. 

In this work, we will investigate the symmetric spin- 
1/2 Ising-Heisenberg diamond chain, which extends the 
class of lattice-statistical models considered in Refs. l32l 
35) by including the Ising four-spin interaction. Beside the 
ground-state analysis, our main goal is to particularly ex- 
amine the effect of the four-spin interaction on the mag- 
netization scenario as well as the thermodynamics of the 
Ising-Heisenberg diamond chain. Note that the considered 
multispin interaction does not violate the standard com- 
mutation relations for cluster Hamiltonians of two differ- 
ent diamond-shaped units, and hence all the results for the 
investigated extended version of the Ising-Heisenberg di- 
amond chain will be attained by combining the general- 
ized decoration-iteration mapping transformation 14611471 
I48H491 and the transfer-matrix technique 144U45I1 . 

The organization of this paper is as follows. In Sec- 
tion |2j we will define the symmetric spin- 1/2 Ising- 
Heisenberg diamond chain with the four-spin exchange 
interaction and then, the most important steps of an ex- 
act analytical treatment will be reviewed. Section [3] deals 
with a discussion of the most interesting numerical results 
for the ground state, the magnetization process, the sus- 



ceptibility and the specific heat. Finally, some concluding 
remarks are drawn in Section[4] 

2 Model and its exact solution Let us consider 
an one-dimensional lattice of A*" inter-connected diamonds 
(see Fig. [TJ defined by the Hamiltonian H = X^feLi ^fc> 
where 



U k = J H A{St k _ x S., 



3k 



°3k- 



1^3k) + S-ik-iS; 



3k 



+ M S 3k-l + S 3k)(^3k-2 + Vtk+l) 

+ K ^3k-1^3k^3k-2^3k+l ~ ^n(^3k-l + ^3k) 



3k-2 



'3k+ 



i)A 



(1) 



Here, Sj ( 7 



j k K j — x,y, z) and a z k denote spatial compo- 
nents of the spin-1/2 operators, the parameter Jn(A) la- 
bels the anisotropic XXZ interaction between the nearest- 
neighbouring Heisenberg spins, in which the parameter A 
allows to control the interaction Jh between the easy-axis 
{A < 1) and easy-plane (A > 1) type, as well as, to ob- 
tain the Ising model as a special limiting case when assum- 
ing A = 0. The parameter J\ denotes the Ising interac- 
tion between the Heisenberg spins and their nearest Ising 
neighbours, and the parameter K describes the Ising four- 
spin interaction between both Heisenberg spins and two 
Ising spins of the diamond-shaped unit. Finally, the last two 
terms determine the magnetostatic Zeeman's energy of the 
Ising and Heisenberg spins placed in an external magnetic 
field H\ and iJn oriented along the z-axis, respectively. 

Exact solution for the considered Ising-Heisenberg di- 
amond chain can be achieved by following the same proce- 
dure as in our recent work I351 . where we have solved ex- 
actly the simplified version of this model (the model given 
by the Hamiltonian (Q]), but without the term describing 
the four-spin interaction). Indeed, by taking into account 
a validity of the commutation relation between the clus- 
ter Hamiltonians of different diamond units, the partition 
function Z of the investigated Ising-Heisenberg diamond 
chain can be partially factorized: 



/V 



-0U k 



(2) 



2 = E%}^ = EI1 T ^ 

{a k } {<j k } k=l 

Here, (5 = l/(k-QT) (fee is Boltzmann's constant, T is 
the absolute temperature) and N is the total number of 
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the Ising spins. The symbols X){ CTfc } an ^ Tr{s fc } denote a 
summation over spin degrees of freedom of all Ising and 
Heisenberg spins, respectively, while Tr^ means the partial 
trace over spin degrees of freedom of two Heisenberg spins 
from fcth diamond unit. It is obvious from the structure of 
the partition function (0 that it is necessary to compute 
the last partial trace Tr k e~P Hk in order to proceed fur- 
ther with a calculation. For this purpose, it is useful to 
pass to the matrix representation of the Hamiltonian (0 
in the orthogonal basis constructed by the eigenstates of 
the operator S'| fc _ 1 5'| fc : {|+, +) 3 fc-i,3fc. |+, ->3fe-i,3fc. 
|-,+)a*-i,3fc. |-, ->3fe-i,3fe}, where the relevant ket 
vectors | + (-), +(-)} 3 fc-i,3fc = I + (-))sk-i| + (-)>3k 
correspond to the appropriate combinations of the spin 
states 5. 



actually arrives at the following expression for the corre- 
sponding partition function 



Z lG (f3,R,H IC ) 



,JV/3i?./4 



cosh(/3ifi C /2) 



sinh(/3ffic/2) +e-P R 



N 



(8) 



At this stage, exact results for all physical quantities, 
which are important for understanding of the magnetic 
behaviour of the investigated spin system, follow straight- 
forwardly, whereas the Helmholtz free energy of the spin- 
1/2 Ising-Heisenberg diamond chain may be connected to 
the Helmholtz free energy of the uniform spin- 1/2 Ising 
chain (Jrc = — fcB^lnZic) through the relation 



3k-l 



±1/2 and Sf fc = ±1/2 of two Heisenberg J 7 = J" IC — Nk B T In A. 



spins at (3k — l)th and 3fcth decorating sites, respectively. 
The diagonalization of the Hamiltonian (0 transcribed in 
the matrix representation yields the four energy eigenval- 
ues: 

^l( cr 3fe-2 7 cr 3fe+l) = J h/4 + #<4fc-2°3fc+l/ 4 - H VL 

+ (Ji-H I /2)(al k _ 2 +al k+1 ), (3) 
£2(<rh-2, °Zk+i) = ^h/4 + Ka z 3k „ 2 al k+1 /4 + H H 

-(^ + #1/2X^-2+^+1); W 
f 3,4(o-3 fc _ 2 ,CT3 fc+1 ) = -J H /4- Ka^ k _ 2 a^ k+1 /4 

±J H A/2-H l (l k _ 2 +al k+1 )/2. (5) 

The eigenenergies (0-(0 can be straightforwardly used 
to obtain the last trace emerging in the expression of the 
partition function (0. The resulting expression immedi- 
ately implies the possibility of performing the general- 
ized decoration-iteration mapping transformation given 
by Eq. (3) in Ref. 051 with the modified function 
G[(3Ji{al k _ 2 + <J* k+1 ),PK<jl k _ 2 a* k+1 }: 

G(x, y) = 2e-P J * /A - y/i cosh(x - f3H H ) 

+ 2e^ J H/ 4 +^ 4 cosh(/3J H ZV2). (6) 

Finally, a direct substitution of the mapping transforma- 
tion given by Eq. (3) in Ref. [35| into the expression 
yields to the simple relation between the partition func- 
tion Z of the investigated mixed-spin Ising-Heisenberg di- 
amond chain and the partition function Zjc of the uniform 
spin- 1/2 Ising linear chain with the nearest-neighbour cou- 
pling R and the effective magnetic field H\c 

Z(f3,J h J }i ,K,A,H h H }i ) = A N Z lc (f3,R,H lc ). (7) 

The mapping parameters A, R and H\q emerging in Eq. (0 
can be obtained from the 'self-consistency' condition of 
the applied decoration-iteration transformation (for more 
details see Ref. J35]). The mapping relationship be- 
tween the partition functions Z and Z\q completes the ex- 
act calculation of the partition function, since the partition 
function of the uniform spin- 1/2 Ising chain can simply 
be calculated within the framework of the transfer-matrix 
method [44,45 1. In the thermodynamic limit N — > 00, one 



(9) 



2.1 Magnetization and correlation functions The 

sublattice magnetization m\ and mn reduced per one Ising 
and Heisenberg spin, respectively, can easily be calculated 
by differentiating Eq. (0 with respect to the particular 
magnetic field acting on Ising and Heisenberg spins, re- 
spectively: 



mi : 



m H 



1 / dT \ _ 1_ 



N 
1 



2N\dH H J T 



10 



N \dH IC ) T 



mic, 



(10) 



7 2 



mic+eic 



j3 



-£IC 



2eic ■ (ID 



In the above, the parameters m\c an d £ic denote the 
single-site magnetization and the correlation function 
between the nearest-neighbour spins of the uniform 
spin- 1/2 Ising chain 144 1145 II and the coefficients Jjj^ 
(i = 1,2,3) mark the expressions Jh k = -jp— lnGi 

with G t = G(PJ Y ,PK/A), G 2 = G(-0J h 0K/4) and 
G 3 = G(0, -0K/4). In view of this notation, the total 
magnetization normalized per one spin of the diamond 
chain can be expressed as 

m = -(mi + 2m H ). (12) 
o 

It is worthwhile to mention that both sublattice magneti- 
zation listed above can also be calculated by combining 
the mapping relation with the exact mapping theo- 
rems developed by Barry et al. [ 62, 63, 64 1 and the gen- 
eralized Callen-Suzuki spin identity [65,66 67 1. This al- 
ternative approach is more general, since it enables to 
evaluate the ensemble average for any combination of 
random spin variables involved in the total Hamiltonian 
of the system. For example, one easily obtains the ex- 
act expressions for other important physical quantities 
such as the pair correlation functions Cff = (<7| fc _ 2 ), 

Qffl = Wlfc-l^lfe)' ^HH = {^3k-1^3k)> an( l Qh = 

(^3k-2^3k-i) as we H as tne four-spin correlation function 
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(^3fe-2^3fc+i%fe-1^3fe) in Ms Wa Y : 



Chh = ^i"/4 + Ki mic + Ki e IC , 
Chh = + K.2 mic + £3 ei C , 

= /C 3 /2 + (/C+ + £3 ) mic + 2/C 3 eic 
QS = /Cf/16 + /Cimic/4 + /C+ei C /4 



(13) 
(14) 
(15) 
(16) 
(17) 



The symbol (...) denotes the standard canonical average 
performed over the ensemble defined on the investigated 
Ising-Heisenberg diamond chain and the coefficients K,f 
and JCi (i = 1. 2, 3) emerging in Eqs. (TT~4T > — (fl~7T> mark the 
expressions JCf = F^Ji^K/4) + Fi(-j3J h /3K/4) ± 
2^(0,-/3^/4) and /Q = F^Ji, pK/4)-Fi(-pj It 0K/4) 
with 



-/3J H /4-j//4 



cosh(a; — /3-Hh) 



^ JH/4+y/4 cosh(/?J H Zy2)]/ [2G(x, y)], 
^j h /4+j,/4 &&(pj n A/2)/ [2G[x, y)\, 

3 -,3JH/4-y/4 sinh ^ _ I l£G(x, y)] . 



Fi(x,y) ■■ 

F 2 (x,y) 
F 3 (x,y) ■■ 

2.2 Thermodynamics Before concluding this sec- 
tion, it is worthy of notice that several basic thermody- 
namic quantities such as the susceptibility \ an d the spe- 
cific heat C can also be readily derived from the Helmholtz 
free energy (0 by using the standard thermodynamic rela- 
tions 

However, the final expressions for the quantities \ an d C 
are too cumbersome, therefore we do not write them here 
explicitly. 

3 Results and discussion In this section, we will 
focus our attention to the most interesting results to be 
obtained for ground state and basic thermodynamic quan- 
tities of the investigated diamond chain. Even though all 
the results derived in the preceding section are general and 
hold regardless of whether ferromagnetic or antiferromag- 
netic exchange interactions are assumed, we will consider 
both the exchange parameters J\ and Jh to be antiferro- 
magnetic ( J\ > 0, Jh > 0) only, since it can be ex- 
pected that the magnetic behaviour of the model with the 
antiferromagnetic interactions in the external longitudinal 
magnetic field should be much more interesting compared 
with its ferromagnetic counterpart. Moreover, our subse- 
quent analysis will be restricted to the special case with 
the uniform magnetic field iJjj = H 1 = H acting on the 
spin-1/2 atoms, which physically corresponds to the situa- 
tion with the equal g-factors for both kinds of the magnetic 
atoms. Under these assumptions, we can set the Ising in- 
teraction Ji as the energy unit and introduce the following 
set of dimensionless parameters: a = Jh/Ji, 0:4 = K/J\, 
h = H/Ji, and t = ksT/Jj, as describing the strength 
of the Heisenberg interaction normalized with respect to 



the Ising interaction, the strength of the four-spin interac- 
tion normalized with respect to the Ising interaction, the 
relative strength of the external magnetic field and the di- 
mensionless temperature, respectively. 

3.1 Ground-state phase diagrams First, let us con- 
struct the ground-state phase diagrams for the antiferro- 
magnetic spin-1/2 Ising-Heisenberg diamond chain with 
the four-spin interaction under the assumption of the fixed 
exchange anisotropy parameter A = 1, which clarify the 
magnetic behaviour of the system without the external 
magnetic field and the magnetic behaviour of the system 
in the external magnetic field. The ground state phase di- 
agram for the spin-1/2 Ising-Heisenberg diamond chain 
with the four-spin interaction without the external mag- 
netic field is illustrated in Fig. |2a). As one can see, three 
different phases FRIi, QFI and QAF appear in the ground 
state of the system without the external magnetic field due 
to a mutual interplay between the exchange interactions 
Jh, Ji, K and the fixed exchange anisotropy A = 1. They 
spin arrangements can be unambiguously characterized as 
follows: 

N 

I FRIi) = II I->3fc-aI+,+>3*-i,a*, 
fe=i 

mi = -1/2, m H = 1/2, Cfi z = 1/4, 

= 1/4, Cgf, = 0, C§ = -1/4, Qfi = 1/16, 
fraii =Jk/4- Ji + K/16 -H/2; (19) 



|QFI) = ni+) 3fc _ 2 i(|+ ! -)-|-, 



fc=i 



3fe-l,3fc 



mi = 1/2, m H = 0, Cg = 1/4, 

= -1/4, qf H - -1/4, Cfi = 0, QS = -1/16, 
£qfi = -Jh/4 - JuA/2 - K/16 - H/2 ; (20) 



N 



|QAF) = niH 



fc=l 



3fe " 2 ^ ^(|+,-)-K 



)3k-2 



3fc-l,3fc 



mi = 0, m H = 0, Cff = -1/4, 

= -1/4, qf H = -1/4, Cfi = 0, QS 
£qaf = - Jh/4 - Jh^/2 + K/16 . 



1/16, 



(21) 



In above, the ket vectors |±) after relevant products de- 
scribe the spin states a z = ±1/2 and S z = ±1/2 of the 
Ising and Heisenberg spins, respectively. The wave func- 
tion of the phase QAF is written with the help of the ex- 
pression (— ) 3fc ~ 2 £ {+,—}, where k — 1,...,N. Obvi- 
ously, all three phases meet at one triple point given by the 
condition 



T = [04, a] = 



0, 



1 



(22) 



As one clearly sees from Eqs. ( fT9b , the phase FRIi rep- 
resents the semi-classically ordered ferrimagnetic phase 
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in other two phases QFI and QAF one observes the effect 
of quantum fluctuations on Heisenberg bonds. Namely, the 
pairs of decorating Heisenberg spins reside at a quantum 
superposition of spin states described by the antisymmetric 
wave function (| + , — ) — | — , +)) 3 fe_i j,J ^ m ^ om me 
phases QFI and QAF. In the former phase QFI, all nodal 
Ising spins occupy the spin state a z = 1/2, while in the 
latter one QAF, one finds a perfect antiferromagnetic or- 
der in the Ising sublattice. Moreover, one may conclude 
from the location of QAF in phase diagrams shown in 
Figs. |2ja) and (b) that this phase is a result of a mutual 
interplay between the quantum fluctuations arising from 
the antiferromagnetic Heisenberg interaction Jh and the 
ferromagnetic four-spin interaction K. Actually, the phase 
QAF appears in the ground state only when a > 1 and 
Q!4 < 0, regardless of whether the external magnetic field 
h is zero or non-zero. Furthermore, the phase diagram 
plotted in Fig. [2jc) suggests that in response to the applied 
magnetic field h there also may arise the ferrimagnetic 
phase FRI2 and the saturated paramagnetic phase SPP in 
the ground state besides the aforementioned FRIi , QFI 
and QAF phases: 



JV 



|FRI 2 > = J] l(-) 3 ^ 2 >3 fe -2| + ,+) 3fc -l,3 fe , 
k=l 

nil =0,rriH = l/2,Cff = 

c& = 1/4, c£h = 0, eg 

£fri 2 = Jh/4 - K/16 - H 



-1/4, 
0, Qm 



1/16, (23) 



3 



a =0.5 




4 SPP 




FRI 2 




QFI 


FRIj 





(c) a 

Figure 2 Ground-state phase diagrams of the antiferro- 
magnetic spin- 1/2 Ising-Heisenberg diamond chain with 
the four-spin interaction (a) for the system without the ex- 
ternal magnetic field and (b),(c) for the system in the exter- 
nal magnetic field with the fixed 04 = —0.5, 0.5, respec- 
tively, when A = 1. 



with the antiparallel alignment between the nearest- 
neighbouring Ising and Heisenberg spins. By contrast, 



/V 

|SPP> = [] l+>3fc-a|+ J +) 3 fc-i,3fc, 
fe=i 

mi = 1/2, m H = 1/2, C{{ = 1/4, 

C^h = 1/4, = 0, C£ = 1/4, Qg = 1/16, (24) 

£ SPP = J H /4 + Ji + K/16 - 31T/2 . 

It can easily be understood from Eqs. d23l that in the FRI2 
state, the pairs of Heisenberg spins occupy the ferromag- 
netic state |+, +)3fe-i 3fc, while the nodal Ising spins are 
ordered antiferromagnetically against each other. It is clear 
that this phase appears as a result of a mutual interplay be- 
tween applied magnetic field h and the antiferromagnetic 
four-spin interaction K, since it is stable in the ground state 
only when h > 1 and 0:4 > [see Fig. [2c)]. Finally, as one 
can expect, the system always undergoes a phase transition 
towards the fully saturated paramagnetic phase SPP, where 
all Ising and Heisenberg spins are completely aligned to- 
wards the external-field direction [see Eqs. (l24l l in high 
magnetic fields. 

3.2 Magnetization process Next, let us turn our at- 
tention to the magnetization process of the investigated 
diamond chain. For this purpose, the total and sublat- 
tice magnetization against the external magnetic field are 
plotted in Fig. [3] for several values of the interaction pa- 
rameters a and 04. In agreement with the ground-state 
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Figure 3 The zero-temperature total and sublattice magnetization against the external magnetic field for the fixed four-spin 
interaction Q4 = —0.5 (left column), a± = 0.5 (right column) and several values of the interaction parameter a. 
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phase diagrams shown in Figs. [2b) and (c), the plotted 
zero-temperature magnetization curves reflect altogether 
six different sequences of field-induced phase transitions: 
FRIi-SPP, QAF-FRIi-SPP, QAF-QFI-SPP, FRIi-FRI 2 - 
SPP, QFI-FRI2-SPP and QFI-SPP. Moreover, it is clear 
from Fig. [3]that one and/or two intermediate plateaus can 
be found in magnetization curves before the total magne- 
tization m tends towards its saturation value. It is worth- 
while to remark that the identified magnetization plateaus, 
appearing at 0, 1/3 and 2/3 of the saturation magnetization 
when normalizing the total magnetization with respect 
to its saturation value, satisfy the Oshikawa-Yamanaka- 
Affeck rule p(S u — m) £ integer l68l . which has been 
proposed as a necessary condition for the formation of 
quantized plateaus (p is a period of the ground state, S u 
and m are the total spin and total magnetization of the 
elementary unit). 

3.3 Thermodynamics Next, let us turn our atten- 
tion to temperature dependencies of some thermodynamic 
quantities. Figure |4] depicts the magnetic susceptibility 
multiplied by the temperature (xt) as a function of the 
temperature in the zero field for several values of the four 
spin interaction 04 when a = 0.9 and 1.2. As can be 
seen, the thermal dependencies of the product \t either 
look like those that are typical for ID quantum ferrimag- 
nets [69 70 71 1 or have an antiferromagnetic character. 
Actually, when the values of the four spin interaction 04 
and the interaction parameter a are chosen so that the 
phases FRIi or QFI constitute the ground state, then the 
quantity exhibits a round minimum upon cooling and then 
exponentially diverges under further temperature suppres- 
sion [see for instance the curves for Q4 = —0.5, 0.1, 0.9 in 
Fig-ffia) and for 0:4 = 0.1, 0.5 in Fig.|5b)]. These depen- 
dencies clearly reveal a ferromagnetic-antiferromagnetic 
crossover. The temperature variation of xt is generally a 
monotonically decreasing function for ferromagnets and 
monotonically increasing function for antiferromagnets 
as the temperature increases. In this respect, the marked 
low-temperature divergence of xt emerges because of the 
ferromagnetic excitations of spins from the magnetically 
ordered ground state, while the monotonous increase of 
Xt observed at higher temperatures is the result of the 
temperature-induced antiferromagnetic spin excitations. 
On the other hand, the quantity xt always exponentially 
tends to zero in the limit t — s- when the interaction param- 
eters «4 and a are chosen so that the phase QAF constitutes 
the ground state [see the curves for 04 = —0.1 and —0.5 
in Fig. Sfb)]. Finally, other notable thermal variations of 
the magnetic susceptibility multiplied by the temperature 
occur for the boundary values of the four-spin interaction 
«4 = 0.8 and 0.0 that correspond to the phase boundaries 
FRIi-QFI and QAF-QFI in the ground state, respectively. 
Here, the magnetic susceptibility diverges as t^ 1 at low 
temperatures and the product xt tends to the value 1/4 
in the zero-temperature limit, which can be explained in 
terms of the Curie law. 
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Figure 4 The magnetic susceptibility multiplied by the 
temperature as a function of the temperature in the zero 
magnetic field for several values of the four-spin interac- 
tion CK4 when (a) a = 0.9 and (b) a = 1.2. 



Another quantity, which is important for understand- 
ing of thermodynamics, is the specific heat. Some temper- 
ature variations of the zero-field specific heat are plotted 
in Fig. [5] for the fixed a = 1.2 and several values of the 
four-spin interaction 04. Note that values of the exchange 
parameters a and 04 are chosen so as to match the re- 
gion near the triple point T given by Eq. (l22l . where the 
phases FRIi, QAF and QFI coexist in the ground state [see 
Fig. [2a)]. As one sees from Figs. [2 a) and (c), if the four- 
spin interaction is chosen to be sufficiently far from the 
triple point, then temperature dependencies of the specific 
heat with one or two round maxima are observable. Obvi- 
ously, besides the usual Schottky-type maximum detected 
at high-temperature part of the specific heat, there also ap- 
pears an additional second maximum located in the low- 
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Figure 5 Temperature variations of the zero-field specific heat for a = 1.2 and several values of the four-spin interaction 
0:4. The right column illustrates in detail the regions rather close to the triple point, where the phases FRIi, QFI and QAF 
coexist in the ground state. 



temperature part of these curves [see the curves labeled as 
ct4 = —1.0 in Fig. [5J a ) an d as a i = 2.0 in Fig. |3c)]. 
As expected, the low-temperature peak becomes more pro- 
nounced, the closer value of the four-spin interaction «4 is 
selected to the triple point [see e.g. the curves labeled as 
ay = 2.0 and 0.5 in Fig.[5jc)]. These observations suggest 
that the double -peak structure in the specific heat curves 
originates from thermal excitations between the ground- 
state spin configuration and the ones close enough in en- 
ergy to the ground state. Apart from these rather trivial 
findings, a remarkable triple-peak specific heat curve can 
also be found when the region rather close to the triple 
point is considered [see the curves labeled as ay = —0.1 
in Fig. [2a) and ay = 0.1 in Fig. [2c)]. Besides two afore- 
described peaks, there also occurs an additional third peak 
to be located at very low temperatures. There are strong in- 
dications that an appearance of this sharp maximum is the 
result of strong thermal excitations between the ground- 
state configuration of Ising spins in QAF and the ones in 
QFI and FRIi. In accordance with this statement, the pecu- 
liar third maximum gradually shifts towards higher temper- 



atures with retiring from the triple point until it coalesces 
with the second low-temperature peak [see Figs. [5jb) and 
[2d)]. Finally, it is worth noting that remarkable thermal 
dependencies of the specific heat with triple-peak struc- 
ture can also be detected in the region sufficiently close 
to the triple point T when very small external fields are ap- 
plied to the system. The similar triple-peak specific heat 
curves have been recently observed in the simpler version 
of this model without the four-spin interaction when ap- 
plying small external field to the system driven by the spin 
frustration into the disordered ground state ll35ll . 



4 Concluding remarks In the present paper, the 
ground-state properties, magnetization process and ther- 
modynamics of the symmetric antiferromagnetic spin- 1/2 
Ising-Heisenberg diamond chain with the Ising four-spin 
interaction has been particularly investigated within the 
framework of the generalized decoration-iteration map- 
ping transformation. We have shown that the investigated 
diamond chain with the Ising four-spin interaction has 
five different ground states FRIi, QFI, QAF, FRI2 and 
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SPP, providing the antiferromagnetic pair interactions (the 
Ising as well as the Heisenberg ones) and the fixed ex- 
change anisotropy A = 1. We have also found a rigorous 
evidence for one and/or two quantized plateaus in mag- 
netization curves. Our further goal was to shed light on 
thermodynamics of the studied system. The results ob- 
tained for the magnetic susceptibility clearly demonstrate 
that the thermal dependencies of the zero-field suscepti- 
bility multiplied by the temperature either look like those 
that are typical for ID quantum ferrimagnets or have an 
antiferromagnetic character. The most interesting result to 
emerge in the present study is the triple-peak structure in 
the specific heat curves that appears when assuming the 
ground-state region rather close to the triple point, where 
the phases FRIi, QFI and QAF coexist. 

Last but not the least, it should be mentioned that even 
though our theoretical investigation of the symmetric spin- 
1/2 Ising-Heisenberg diamond chain with the four-spin in- 
teraction has been mainly aimed at providing a deeper in- 
sight into the ground-state properties, magnetization pro- 
cess and thermal dependencies of the susceptibility and 
specific heat, it could be quite interesting to explore also 
an enhanced magnetocaloric effect in this diamond chain. 
Beside this, other particular case of this model with the fer- 
romagnetic Heisenberg interaction ( Jh < 0) deserves the 
attention, since there are some fundamental differences be- 
tween magnetic behaviour of the hybrid Ising-Heisenberg 
models with distinct nature of the Heisenberg interaction 
(see Refs. ||29ll54ll55ll58ll59l). Finally, several further in- 
teresting extensions of the present version of the spin- 1/2 
Ising-Heisenberg diamond chain come also into question. 
For instance, it is possible to extend the present model 
by including the cyclic Heisenberg four-spin interaction 
between the Ising and Heisenberg spins of the diamond- 
shaped unit. In this direction will continue our next work. 
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